Photonics meets excitonics: natural and artificial molecular aggregates 
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Organic molecules store the energy of absorbed light in the form of charge-neutral molecular 
excitations - Frenkel excitons. Usually, in amorphous organic materials, excitons are viewed as 
quasiparticles, localized on single molecules, which diffuse randomly through the structure. How- 
ever, the picture of incoherent hopping is not applicable to some classes of molecular aggregates - 
assemblies of molecules that have strong near field interaction between electronic excitations in the 
individual subunits. Molecular aggregates can be found in nature, in photosynthetic complexes of 
plants and bacteria, and they can also be produced artificially in various forms including quasi-one 
dimensional chains, two-dimensional films, tubes, etc. In these structures light is absorbed collec- 
tively by many molecules and the following dynamics of molecular excitation possesses coherent 
properties. This energy transfer mechanism, mediated by the coherent exciton dynamics, resembles 
the propagation of electromagnetic waves through a structured medium on the nanometer scale. 
The absorbed energy can be transferred resonantly over distances of hundreds of nanometers before 
exciton relaxation occurs. Furthermore, the spatial and energetic landscape of molecular aggregates 
can enable the funneling of the exciton energy to a small number of molecules either within or 
outside the aggregate. In this review we establish a bridge between the fields of photonics and 
excitonics by describing the present understanding of exciton dynamics in molecular aggregates. 
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I. INTRODUCTION 

Advances in nanotechnology supported by our under- 
standing of material properties on the microscopic level 
persistently drive the field of photonics to the nanome- 
ter scale. With the development of photonic crystals^ 
and semiconductor optical cavities 3,4 the size of optical 
devices, usually limited by light diffraction in free space, 
can been scaled down to hundreds of nanometers just by 
exploiting the material's dielectric constant. Advances 
in the design of plasmonic metamaterials^— permitted 
this limit to be pushed even further down to the order 
of tens of nanometers. The natural question which arises 
is: what will be the next limit and how can it be ap- 
proached? 

In this review we consider one of possible pathways 
towards reaching the goal of optical devices on the 
nanometer scale: using molecular aggregates as photon 
processing elements on the true nanometer scale, see 
Fig. [TJ Molecular aggregates are assemblies of molecules 
held in place by non-covalent interactions. These sin- 
gle molecules mostly keep their nuclear and electronic 
structure. Examples are molecular crystals as well as 
nanoscale self- assembled structures, molecular films and 
light harvesting systems in photosynthesis. If the lowest 
allowed electronic transitions in the composing molecules 
are within the visible part of the optical spectrum and 
if the molecules have large absorption and fluorescence 
cross-sections, the interaction between molecular elec- 
tronic transitions is strong enough to transfer the light 
absorbed from one molecule to the others, via a resonant 
de-excitation/excitation process. 

This mechanism is closely related to the near-field en- 



ergy transfer between classical dipoles (antennas^—. 
In this sense, the exciton dynamics in molecular aggre- 
gates resembles the propagation of light in a metamate- 
rial where dye molecules play the role of functional ele- 
ments. In contrast to excitons in inorganic semiconduc- 
tors, where free charge carrier mobility determines the 
exciton transport, in aggregates, electrons remain local- 
ized on molecules while the excitation is transferred. 

A canonical example of molecular aggregates that pos- 
sess coherent exciton properties are J-aggregates. These 
are aggregates of fluorescent molecules, discovered about 
80 years ago independently by Scheib o 14 i 15 and Jelle)*i£. 
They observed that the formation of aggregates is ac- 
companied by drastic changes in the optical properties. 
The broad absorption line corresponding to the molec- 
ular excitation is shifted to the red side of the spec- 
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FIG. 1: Schematic illustration of the length scales character- 
izing the different physical phenomena present the different 
fields discussed in this review. 



trum and becomes much narrower. Later, these types 
of molecular aggregates were named J-aggregates after 
Jelley. J-aggregates represent only one class of molecu- 
lar aggregates. For instance, there exist aggregates where 
the absorption line is blue-shifted, so called H-aggregates. 
Some aggregates exhibit several J-bands 17 and some have 
both J-and H-bands^ 8 . In aggregates that consist of non- 
equivalent monomers even more complicated absorption 
structures can be found 19 . 

Molecular aggregates also appear as functional units in 
nature. For instance, they form the absorbing and energy 
transferring parts of light-harvesting complexes in plants 
and some types of bacteria and algae—. In these com- 
plexes the exciton absorbed by the antenna aggregate has 
to be funnelled to the reaction center - the part where 
the exciton energy is used to create free charges to be 
employed in a chemical reaction. The efficiency and ro- 
bustness of light absorption and exciton transfer in light 
harvesting complexes may be crucial for the survival of 
photosynthetic organisms under evolutionary pressure. 

Shortly after the discovery of the self- assembled or- 
ganic dye aggregates the close connection of these arti- 
ficial structures to the natural light harvesting systems 
was recognize d 20 ! 21 and the exciton model of Frenkel 9 -* 2 ^ 
has been used to explain the observed changes in op- 
tical properties and the transfer of energy along the 
aggregate 21 . 

In Frenkel's exciton theory, which is based on the 
classical resonance interaction theory of Holtsmark-, the 
electronic excitation in the aggregate is not confined to 
a single monomer, but it is coherently delocalized over 
many monomers in the form of "excitation waves" . From 
superpositions of these excitation waves "excitation pack- 
ets" can be formed, which describe the coherent motion 
of (localized) excitations 22 . Already in these early works 
it was established that coupling to internal and external 
vibrational modes, imperfections of the aggregates and 
disorder induced by the environment strongly influences 
the "coherence size" of the exciton waves and modifies 
the optical and transport properties. 

When the interaction with the environment (internal 
vibrations, solvent, etc.) is much stronger than the 
resonant excitation transfer interaction, the excitation 
becomes more or less localized on one monomer and 
the transfer is no longer described by coherent exciton 
motion but it becomes an incoherent hopping process. 
Forster derived an elegant formula for the rate constants 
for transport of excitation from one monomer to the 
othe r 23 ! 24 . This rate is proportional to the overlap of 
the donor emission spectrum and the absorption spec- 
trum of the acceptor molecule and depends on the inverse 
sixth power of the distance between donor and acceptor. 
Typically, in molecular aggregates, one is in a regime 
in-between these two extreme cases, i.e. the transport is 
neither fully coherent nor incoherent. This complicates 
the theoretical modeling and the interpretation of exper- 
iments. 

The remarkable optical and transport properties of 



molecular aggregates have led to a variety of appli- 
cations. Right from the beginning, molecular aggre- 
gates were employed as wavelength selective sensitiz- 
ers in photograph y 25 ! 26 . Recent applications include 
the measurement of membrane potentials 27 or the de- 
sign of colorants 2 ^. Some molecular aggregates form 
self- assembled supramolecular flexible fluorescent fibers 29 
which may have applications in thin-film optical and 
optoelectronic devices, for example by employing opti- 
cal bistability of J-aggregates^ -. Molecular aggregates 
can be also utilized for sensing applications. For in- 
stance, it has been demonstrated that a large absorption 
cross section combined with fast exciton diffusion may 
be used to enhance fluorescence from a small number 
of dye molecules adsorbed or embedded in an aggregate 
fil m 31 ! 32 . Due to fast exciton diffusion within the aggre- 
gate the excitons explore the aggregate. Once they find 
a molecule with an electronic transition close to the J- 
band the exciton can be transferred inelastically to the 
adsorbed molecule. Moreover, dye aggregates might play 
an important role in the development of efficient, low-cost 
artificial light harvesting units (like organic solar cells)^ 3 .. 
This review gives a brief, yet by no means complete, 
overview of our experimental and theoretical understand- 
ing of excitons in molecular aggregates. Most of the as- 
pects in this work are covered at the advanced intro- 
ductory level, and for a deeper understanding we refer 
to more detailed studies such a o 19 i 34 ~ — . The review is 
structured as follows: In Sec. [Ill we describe basic ex- 
citonic properties of molecular aggregates starting with 
single molecules and molecular dimers and then introduc- 
ing several examples of artificial and natural molecular 
aggregates. In Sec. 1111} we overview the main experimen- 
tal techniques that allow for probing the structural prop- 
erties of molecular aggregates and excitation dynamics. 
Section [TV] introduces the main theoretical approaches 
utilized for modeling of excitons in aggregates. Section 
El shows how molecular aggregates can be combined with 
photonic and plasmonic structures. Finally, we conclude 
the review in Section EH 



II. 



BASIC PROPERTIES OF MOLECULAR 
AGGREGATES 

A. Properties of monomers 



Molecules in the aggregates largely retain their elec- 
tronic and nuclear structure. Thus, it is natural to start 
our discussion with individual molecules interacting with 
light. The interaction of a single molecule with light is 
schematically illustrated in Fig. [2j If the molecule has an 
optically allowed electronic transitions within the photon 
spectrum it can be excited by absorbing light at the cor- 
responding frequency. The transition time between the 
ground and excited state of the molecule is of the order of 
tens of attoseconds - the time that a single photon inter- 
acts with a molecule. During this short time interval the 
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FIG. 2: Excitation/de-excitation processes in a single 
molecule illustrated on a two-level molecular energy diagram. 
The parabolic surfaces correspond to the ground state \g) and 
the first excited |e) electronic states of the molecule. The hor- 
izontal axis indicates the displacements of nuclei from their 
equilibrium positions. The molecule, initially in the ground 
electronic and vibrational states, is excited by absorbing light 

- blue arrow. During the absorption process the positions of 
the nuclei are not changed. The light absorption also induces 
molecular vibrations. Due to the interaction with the en- 
vironment the molecular vibrations are equilibrated and the 
molecule relaxes to the bottom of the excited electronic state 

- orange wavy arrow. Finally, the molecule relaxes down to 
the ground electronic state by emitting a photon - red arrow 
(fluorescence) , or without photon emission - green dashed ar- 
row (non-radiative process). In the ground state the induced 
vibrations are also equilibrated due to interaction with the 
environment. 



positions of the nuclei in the molecule are not changed. 
Initially the molecule is in its ground state geometry, i.e. 
in a minimum of the electronic ground state potential \g) 
indicated by the blue wave-packet in Fig [2j The equi- 
librium positions of atoms for the excited states, usu- 
ally, are different from those of the ground state. There- 
fore, after the transition to the state |e) the molecule is 
in a transient non-equilibrium state where molecular vi- 
brational modes are also excited (green wave-packet in 
Fig [2]). Then, due to the interaction of the molecule 
with its environment the molecule relaxes towards the 
energetic minimum of the excited state, i.e. the excited 
vibrational modes are relaxed (wavy arrow and orange 
wavepacket in Fig [2j). This relaxation of the molecu- 
lar geometry can also be accompanied by a rearrange- 
ment of the environmental molecules to a configuration 
with lower total energy. At ambient conditions this pro- 
cess occurs over a timescale of hundreds of femtoseconds 
to several picoseconds. Finally, on timescales of tens of 
picoseconds to nanoseconds the molecule relaxes to its 
electronic ground state either emitting a photon (fluo- 
rescence) or transferring energy to other degrees of free- 
dom, for instance vibrations (non-radiative relaxation^). 



Usually, the absorption and emission spectra exhibit a 
progression of peaks stemming from the coupling to vi- 
brational modes with high energy (~150 meV). These 
peaks are broadened by the same order of magnitude, 
due to coupling to a multitude of low energy modes of 
the molecule and the surroundings. Emission takes place 
from the thermally relaxed excited state, which is typ- 
ically located in the low energy wing of the absorption 
spectrum. The relaxation energy is related to the en- 
ergy difference between the maxima of the absorption 
and emission spectra, the so-called Stokes shift. 

The absorption efficiency of a particular electronic 
transition can often be characterized by the correspond- 
ing transition dipole - a matrix element of a dipole opera- 
tor between the ground and the excited molecular states 



d= (e\q-r\g). 



(1) 



Here and in the following we use an arrow symbol to indi- 
cate three-dimensional vectors in real space. The dipole 
operator is characterized by its charge q and the position 
operator f. The transition dipole is not a measurable 
value and is defined up to a complex phase factor (if 
magnetic interactions can be neglected, which is often 
the case, then the wave function can be chosen real and 
d becomes a real vector). The absorption and emission 
strengths of the respective transition are proportional to 
\d\ 2 . In the absorption spectra the transition frequency is 
determined by the energy difference e between the ground 
and excited state. 

Structures and optical spectra of representative 
molecules are shown in Figs. [3] and 2J In Fig. [3] we 
show a TDBC molecule, which can be used as an illustra- 
tive example for a broader class of cyanine dyes. Beside 
the monomeric absorption and fluorescence spectra (in 
solution), spectra of aggregates formed on a glass sub- 
strate are also shown. The lowest electronic transition 
in TDBC is characterized by a large transition dipole 
of the order of \d\ « 10 Debye (Debye units are com- 
monly used for molecular dipoles, 1 Debye « 0.208 e- A 
where e denotes the electron charge), which is aligned 
with the backbone of the molecule. The higher excited 
states in TDBC are well separated from the lowest one, 
which indicates that, for low intensity optical absorp- 
tion and exciton dynamics it is sufficient to take only 
two electronic states into account. The absorption spec- 
trum of the molecule in a solution shows a broad line 
in the green part of the spectrum with a partially re- 
solved vibronic structure. Note also that the emission 
from a monomer is roughly a mirror-image of the ab- 
sorption, where the maximum is shifted to lower ener- 
gies. This is the Stokes shift, which is caused by the 
relaxation/reorganisation after an electronic transition. 
The mirror image of absorption and fluorescence indi- 
cates that the vibrational frequencies and their coupling 
to electronic transitions are very similar in the electronic 
ground and excited state. Upon aggregation the absorp- 
tion and fluorescence lines become much narrower with 
no structure and are shifted to the orange-red color. The 
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FIG. 3: TDBC fluorescent dye that forms ID and 2D molec- 
ular aggregates, (a) Molecular structure, grey and white 
spheres represent carbons and hydrogens respectively, (b) 
Normalized absorption and fluorescence spectra of TDBC 
molecules taken in a solution (solid lines) and in a 2D ag- 
gregated form (dashed lines). 
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FIG. 4: (a) Bacteriochlorophyll (Bchl) molecule - light ab- 
sorbing pigment that is a basic structural element of light- 
harvesting complexes in phototrophic bacteria. Its struc- 
ture is similar to chlorophyll - the photosynthetic pigment 
in plants.; (b) Absorption and fluorescent spectra of Bchl 
molecules. 



Stokes shift in the aggregated form is also much smaller, 
which indicates a reduced coupling to the environment. 
Aggregates composed of TDBC have, for example, been 
produced in solution^ and on surfaces 40 . By changing 
the side groups different geometrical arrangements can 
be achieved^. 

Another representative example of molecules forming 
aggregates is the bacteriochlorophyll (Bchl)- a pigment 
molecule, which is a functional element of photosynthetic 
systems in phototrophic bacteria, see Fig. |4^a). The 
structure of Bchl is similar to that of chlorophyll - the 
photosynthetic pigment in plants. Both are derivatives 
of porphyrin, compelexed with Mg 2+ . The lowest elec- 
tronic transition, which is involved in the energy transfer, 
is the so called Q^-band. The corresponding transition 
dipole lies in the plane of the porphyrin ring and is about 
\d\ ~ 5 Debye. The transition dipole corresponding to 
the next electronic state is denoted as Q x . It is perpen- 
dicular to the Q y dipole, and the associated transition 
strength is much smaller. The Soret band which lies in 
the ultraviolet has stronger absorption but is usually not 
involved in the energy transfer. As in the case of TDBC 
the emission occurs from the lowest excited state and 



the emission spectrum is roughly a mirror image of the 
absorption band corresponding to this state. 



B. The aggregate 

Organic dye molecules can self-aggregate in different 
types of structures. Sometimes the aggregation is driven 
by electrostatic forces pushing molecules to adsorb on a 
surface; sometimes hydrophobic parts of the molecules 
repel the water and tend to collect together such as in 
synthetic tube aggregates^. 

A specific property of molecular aggregates with strong 
inter-monomer couplings is that the absorption and flu- 
orescent spectra substantially differ from the spectra of 
the molecules which form the aggregates while the elec- 
tronic structure of the molecules is not modified. The 
intermolecular distance in an aggregate is large enough 
that the electron tunneling between different molecules 
can be neglected. Therefore, only the Coulomb interac- 
tion between the electronic transitions in monomers is 
responsible for the spectral modifications and the inter- 



molecular energy transfer. This interaction is similar to 
the non-radiative near-field coupling between plasmonic 
structures, which is mediated by virtual photons. It turns 
out that the major contribution of the Coulomb interac- 
tion in the spectra of aggregates is from the excitation 
transfer via resonant transitions. We will refer to this 
main interaction term between the monomers as Forster 
coupling. The Forster coupling, to a good approxima- 
tion, can be described by restricting to dipole transitions 
only (see below Eq. [2]). The interaction is visualized in 
Fig. [5] for a situation in which initially monomer 1 is 
electronically excited and monomer 2 is in the electronic 
ground state. Then, in a resonant process molecule 1 is 
deexcited and simultaneously the second monomer is ex- 
cited. While this time-ordered discrete picture is easy for 
visualisation in reality the exciton energy is transferred 
coherently. For large distances between the monomers 



identical monomers, ignoring nuclear degrees of free- 
dom for the moment. The eigenstates of the uncou- 
pled monomers can be taken as a basis. These states 
are \g)\g), \e)\g), \g)\e), |e)|e). While the states \g)\g) 
and \e)\e) are still suitable to describe the ground state 
and the doubly excited state of the dimer, respectively, 
the two degenerate single exciton states \e)\g) and \g)\e) 
are no longer eigenstates of the coupled system, because 
of the resonant transfer interaction. The correspond- 
ing Hamiltonian is written as H = £i|l)(l| + £2|2)(2| + 
V12 ( 1 1) (2 1 + |2)(1|) (we choose the energy of the monomer 
ground state as zero), where we have introduced the 
shorthand notation |1) = \e)\g) and |2) = \g)\e) or equiv- 
alently in a matrix representation 
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FIG. 5: Illustration of the excitation transfer due to resonant 
near-field interaction between the state \e)\g) and \g)\e). 



the interaction between monomers 1 and 2 can be ap- 
proximated by the dipole-dipole term 

V 12 = ^-(di • d 2 - 3(di • R 12 )(d 2 • £ 12 )), (2) 

where R\ 2 denotes the distance between the monomers, 
R\ 2 is the corresponding direction vector, and the vectors 
d n are the transition dipoles introduced in Eq. (pQ). 

If the intermolecular distance is comparable with the 
size of the molecules, then the interaction between the 
two monomers can no longer be described adequately us- 
ing the point dipole-dipole interaction of Eq. (|2]). Often, 
it is then sufficient to replace the point-dipoles in Eq. (J2j) 
by extended dipoles^. For even higher accuracy more 
elaborate schemes have been developed (see e.g4£). 



For identical energies of excited states s\ = e 2 = e 
the eigenstates of the dimer are superpositions |=b) = 
-W|e)|g) ± |e)|g)) where the electronic excitation is co- 
herently delocalized over both monomers. The corre- 
sponding eigenenergies are e± = e ± V± 2 . Note that the 
magnitude and sign of the interaction depends sensitively 
on the distance and orientation of the two monomers. 
Fig. [6] illustrates how the relative orientations of molec- 
ular transition dipoles change the frequency of electronic 
excitations in a dimer. If the transition frequencies of the 
involved monomers are equal and the transition dipoles 
are parallel, only the transition between the ground and 
the |+) state is optically allowed. It is detuned by the 
energy V\ 2 from that of the non- interacting monomers. 
From Eq. (|2j) one can see that for case (a), where the 
transition dipoles are orthogonal to the distance vector 
between the monomers one has positive V\ 2 ~ \d\ 2 / R\ 2 - 
For the case (b) the interaction is negative with V\ 2 ~ 
— 2\d\ 2 / R\ 2 - Therefore, the alinement of the molecules 
shown in Fig. [6] (a) and (b) will be seen as a blue and red 
shift of the absorption line, respectively. As it has been 
noted before, the states with only one excited monomer 
such as \g)\e) and \e)\g) are not eigenstates of the sys- 
tem anymore. Thus, if one of these state is populated 
initially, the energy will oscillate between them back and 
forth coherently. 

Let us briefly mention that the dimer often appears 
as a first step of aggregation^ and thus has been inves- 
tigated intensively (see e.g4^— ) . This interest in the 
dimer is also due to the fact that many of the more re- 
alistic models for aggregates where the environment and 
the internal vibrations are included can only be solved 
efficiently for the dimer. 



1. The dimer 

Before discussing general molecular aggregates, let us 
describe some basic results for the case of two coupled 



2. Arbitrary aggregates 

The theoretical description of the dimer can easily be 
extended to arbitrarily arranged monomers. As in the 
case of the dimer the total ground state is taken as a 



V n =\d\ z /Rt 2 
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FIG. 6: Interaction between electronic transitions dipoles 
(blue arrows) of two molecules. The directions of the ar- 
rows reflect the relative phase of the transition dipoles. For 
the considered molecular alignments only one of the optical 
transitions is allowed for the dimer (green arrows), (a) The 
optically allowed transition of the dimer is shifted to the blue 
part of the spectrum, which is similar to H- aggregation; (b) 
the optically allowed transition is shifted to the red, which 
correspond to J- aggregation. 



product of states where all monomers are in the ground 
state | g^f g ) = \g)i" -\g)N- We are interested in the 
properties of aggregates in the linear absorption regime. 
Therefore, it is sufficient to take into account only the 
states with at most one electronic excitation on the ag- 
gregate, such as 



\ n ) = \g)i'"\e)n-"\g) 



N, 



(4) 



in which monomer n is electronically excited and all other 
monomers are in their electronic ground state. In this 
one-exciton manifold approximation the Hamiltonian can 
be written as 



TV 



H e = ^2e n \n)(n\ + ^2v nm \n)(m\ 



(5) 



n=l 



where V nm is the Coulomb dipole-dipole interaction, see 
Eq. [2j which transfers excitation from monomer n to m. 

Examples of commonly discussed ID and 2D planar 
structures are shown in Fig. [71 While ID models are 
well studied theoretically (see e.g. Refs^2~— ) and are fre- 
quently used to characterize excitons in self-assembled 
molecular aggregates, it is not clear whether ideal ID 
molecular chains are formed in experiments 34 . Usually, 
they are subsystems of higher dimensional structures 
such as films or crystals. 

The ID model of molecular aggregates is convenient 
because many results can be obtained analyticall y 19 ! 59 . 
For a perfect very long chain with N molecules the eigen- 



states are well described by "exciton waves" 

N 



1 N 



(6) 



n=l 



with the corresponding eigenenergies 



Ej = e + 2V cos 



(¥)■ 



(7) 



In the last equation, for the sake of simplicity, we have 
considered the interactions between nearest neighbors 
only, which are denoted by V. A discussion of finite one- 
dimensional aggregates can be found e.g. in Refs j 19 ! 59 . 

For the perfect linear chain with parallel monomers, 
the largest transition dipole correspond to the excita- 
tion of the state with a minimal number of nodes. In 
the case of J-aggregation, Fig. [7(a), this state is at the 
bottom of the exciton spectrum shifted by about sev- 
eral hundreds of meV from the monomer transition while 
for H-aggregates, Fig. [7(b), this shift is to the blue part 
of the spectrum. In reality, to compare exciton transi- 
tions in aggregates with excitations in single molecules 
one also has to include Van der Waals interaction with 
off-resonant excitation s 19 ! 47 and the vibronic structure 
of molecular excitatio n 19 ! 47 ! 60 . However, the simplified 
model presented above qualitatively describes the exci- 
ton states in aggregates. 

2D molecular aggregates can be formed, for instance, 
as Langmuir-Blodgett films^i, deposited molecule-by- 
molecule on a surfac e - 40 ! 62 , or by spin coating 63 . The 
structure of the aggregate is determined by the molecular 
properties such as their shape, charge, etc., and the as- 
sembly method. For instance, for cyanine dyes^ (TDBC 
is one example of the class of molecules) some commonly 
assumed aggregation structures are brickstone, Fig. [7(c), 
and herringbone, Fig. [7(d). 

3D molecular aggregates with translational symmetry 
are usually known as molecular crystals. Frequently, in 
these structures is possible to identify lower dimensional 
subsystems with a preferred interaction between molec- 
ular transition dipoles 65 . These directions would deter- 
mine anisotropy of exciton transport. While on average 
the exciton transport in a crystal can be incoherent one 
may expect specific directions with large coherent exciton 
derealization. 

The above mentioned structures are only very simple 
examples of molecular aggregates. For instance, when 
driven by hydrophobic and hydrophilic forces molecules 
can aggregate in two-layer tubes of about 10 nanometers 
in diamete r 17 ! 66 or one can use templates like polypep- 
tides or DNA to induce helical structure s 67 ! 68 . 

Because of the huge variety of organic dyes that can ag- 
gregate it is possible to create narrow J-band absorption 
within an arbitrary part of visible and near IR spectrum, 
see Fig [8l 
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FIG. 7: Illustrations of ID and 2D planar molecular aggre- 
gates. Each block corresponds to a monomer forming the 
aggregate, (a) - staircase and (b) - ladder models for ID 
packing; (c) - brickstone and (d) - herringbone 2D packing 
models. 



FIG. 9: Schematic structure of light-harvesting complex in 
green sulfur bacteria. The solar light is absorbed by chloro- 
some antenna aggregates of Bchl molecules. Then, the energy 
in the form of excitons is transferred through the baseplate 
and Fenna-Matthews-Olson (FMO) protein complexes to the 
reaction center, where the charge separation occurs. 
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FIG. 8: Absorption spectra of four different cyanine dyes in 
monomeric and aggregated forms illustrating the variability 
in the aggregate optical properties. Adapted with permission 
from Ref.— . Copyright 2011 American Chemical Society. 



C. Natural aggregates 

In nature, photosynthetic complexes of plants, pho- 
totrophic bacteria and algae are composed of molecular 
aggregates that efficiently transfer excitons. As an ex- 



ample, in Fig. [9] we sketch the energy transfer along the 
light harvesting system of green sulfur bacteria (GSB). 
These GSB can be found e.g. in deep sea. By develop- 
ing a sophisticated light-harvesting antennae structure 
these organisms adapted to survive at very low light 
intensitie o 70 i 71 . Energy absorption and transfer in GSB 
goes through a network of Bchl molecules, Fig. HJ aggre- 
gated in several types of functional structures - chloro- 
some antenna, baseplate, Fenna-Matthews-Olson protein 
complex (FMO), and the reaction center. 

The photons are absorbed by the chlorosome - an or- 
ganelle which is bound to the bacterial membrane and 
has a size of several hundreds of nanometers. The chloro- 
some contains a disordered array of cylindrical or ellip- 
soidal molecular aggregates - antennas. While differ- 
ent types of geometrical and structural disorders compli- 
cate the molecular structure characterization of natural 
aggregates, recent NMR analysis of mutant chlorosome 
antennas 7 ^ suggested that Bchl molecules in them are 
arranged in an array of concentric helical structures il- 
lustrated in Fig. [TOta). The Forster coupling between 
nearest-neighbor Bchl molecules in the chlorosome an- 
tenna complexes is of the order of 100 meV and the pro- 
posed molecular arrangement of the mutant should result 
in the formation of a J-band. Both experimental 7 -^ and 
theoretical studies 7 ^, show that the exciton spreads over 
a single antennae on the timescale of hundreds of fem- 
toseconds. 

Another functional element of GSB light-harvesting 
structure, is the Fenna-Matthews-Olson (FMO) complex, 
a protein which is depicted in Fig. fTTjTc). FMO plays 
the role of a molecular wire transferring energy from the 
chlorosome to the reaction center. It is a trimer con- 
taining 8 Bchl molecules in each subunit. Unlike in the 
chlorosome antenna, Bchl molecules in FMO are held to- 
gether by a protein cage, and also the Forster coupling 




FIG. 10: Examples of molecular aggregates in the light- 
harvesting complex of green sulfur bacteria, (a) Chlorosome 
antenna complex - a cylindrical or ellipsoidal aggregate of 
Bchl molecules, several concentric aggregates are enclosed; 

(b) A Bchl dimer - a unit block of the antenna complex— ^; 

(c) Fenna-Matthews-Olson (FMO) protein complex - an ex- 
citonic "wire"; (d) organization of Bchls in a monomer unit 
of FMO complex-^ 



between the molecules is several times weaker. However, 
from the multiple experimental and theoretical studies it 
is suggested that the exciton states in FMO are partly 
delocalized 75 77 . 



III. EXPERIMENTAL CHARACTERIZATION 
OF EXCITONIC SYSTEMS 

Much effort is devoted to obtain the microscopic struc- 
ture of the various aggregates (which in turn are needed 
for theoretical models), including the arrangement of the 
monomers, their spacial and energetic disorder, etc. Be- 
side these conformational aspects there is also large inter- 
est in the dynamic properties of excitons such as exciton 
relaxation and dephasing rates, diffusion coefficients, etc. 

X-ray crystallography can be used to obtain lattice 
properties of molecular aggregates, if the latter either 
possess an intrinsically long-range order (like molecular 
crystals or two-dimensional monolayers on substrate^ 9 -) 
or can be crystallized. For instance, crystallization of 
photosynthetic light-harvesting complexes 80 had a large 
impact on the understanding of natural photosynthesis^. 
However, most molecular aggregates cannot be crystal- 
lized. In recent years Cryo- Transmission electron mi- 
croscopy has provided valuable information on the ge- 
ometrical structure of many aggregate s 66 i 81 ~ — . 

The classical way to obtain information on the geome- 
try of the aggregate is by optical spectroscopy (in partic- 



ular absorption, linear dichroism and circular dichroism) 
combined with theoretical modeling. For instance, the 
positions, intensities, polarizations and shapes of the ag- 
gregate absorption bands provide information about the 
strength of intermolecular coupling as well as relative ori- 
entations of the monomer a 19 ! 47 ! 82 . Circular dichroism 
spectra can be used to analyze chirality of the struc- 
tures, for instance, in studies of J-aggregate helices and 
tube o 18 i 84 ~ — . Similarly, the emission spectra can be used 
to obtain information on the aggregate structure. These 
optical experiments are often performed under varying 
environmental conditions like solvent, temperature, con- 
centration of monomers, alignment of the aggregates or 
pressure^—. More detailed information including sep- 
aration of homogeneous and inhomogeneous line widths 
and femtosecond exciton dynamics can be obtained from 
non-linear 2D spectr a 73 ! 91 . Beside optical spectroscopy 
of electronic excitations a multitude of various experi- 
mental techniques is used to characterize aggregates, e.g. 
electroabsorption 92 , Fourier transform infrared (FTIR) 
spectroscopy 9 - 3 -, nuclear magnetic resonance NMR 9 -^ and 
Raman spectroscopy 95 . 

With optical spectroscopy one probes the structure of 
energetic levels in the aggregate and phase relations be- 
tween corresponding electronic transitions, from which 
one can infer the excitation dynamics in the aggregate. 
This might be a viable way for small systems, where the 
exciton is delocalized along the aggregate. However, one 
would like to follow the dynamics of the exciton also in 
real space, that means to measure the time-dependent 
probability to find excitation on a certain monomer. Here 
in particular, diffusion constants are of interest, which 
characterize the spreading of the exciton after the initial 
coherences have died out. Exciton diffusion coefficients 
cannot be measured directly using existing experimental 
techniques. The indirect methods include quenching of 
photoluminescenc o 96 i 97 , photocurrent response 9 -^, tran- 
sient grating 9 - 9 - and exciton-exciton annihilationiSSr-iS 3 -. 
Among the state-of-the-art techniques one could men- 
tion the recently developed coherent nanoscopy^i, which 
would allow for spatial resolution comparable with the 
exciton diffusion length. 



IV. MODELS OF EXCITON DYNAMICS 

While the simple electronic model from Section [II] (see 
in particular Eq. (|5])) already allows us to understand ba- 
sic properties of molecular aggregates (for instance, the 
position of the J-band) it is not sufficient to describe im- 
portant features such as the narrowing of the J-band (in 
particular in contrast to the broad H-band with vibronic 
structure). Also, many transport properties, e.g. depen- 
dence on temperatureiS 5 -, cannot be explained. To this 
end it is necessary to include vibrational modes and the 
influence of the environment into the description. This 
can be done in various ways and there exists a multitude 
of theoretical models in the literature where each is best 
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suited for a particular question and/or situation. For 
example, the Forster rate theory 2 ^ works well when the 
timescale of the transfer is very long compared to deco- 
herence and vibrational relaxation of the excitation. In 
this Section we review the most common models. 



A. Static disorder 

Since usually the experiments are performed on en- 
sembles of aggregates where each aggregate experiences 
a slightly different environment, it is common to average 
over different "configurations" of the aggregate. Most 
often the influence of the environment is treated simply 
as inducing random shifts of the transition frequencies 
of the monomers. For early works see e.g i 54 i 106 i 107 . Of- 
ten, it is assumed that each monomer sees "its own" local 
environment and the shifts induced by this local environ- 
ment are uncorrelated from the shifts of the neighboring 
monomers and disorder of the coupling is ignored. This 
assumption could be called "uncorrelated diagonal disor- 
der" . Sometimes also the positions or couplings of the 
monomers are treated as random variables (off-diagonal 
disorder) 107 i 108 . For each disorder realisation the exciton 
states are typically no longer delocalized over the full ag- 
gregate as in a perfect chain. This localisation becomes 
stronger for larger values of disorde r 54 ! 107 i 109 . In the lit- 
erature, the influence of various forms of correlations has 
been discussed, see e.g j 54 i 11Q i m . To account for tem- 
perature a variant of the open system model described 
below is adaptecUi^^ii^, where the environment of the 
disordered chain can lead to scattering and relaxation 
between the localized exciton states. Such a model has 
been successfully used to describe optical and transfer 
properties of some aggregate a 113 i 115 i 116 . In particular at 
low temperatures the increase of mobility with increasing 
temperature could be explained. 

In the above formulation the influence of internal vi- 
brations of the monomers is neglected. This seems to 
be a good approximation for the J-band where one can 
show that the electronic excitation is to a large extend 
decoupled from vibrational mode a 55 i m i 117 . 



B. Dynamic disorder, Haken-Strobel-Reineker 
model 



Fluctuations of the monomer transition frequencies or 
the intermolecular couplings that are fast compared to 
the exciton transfer time-scale are usually referred to as 
dynamic disorder. Dynamical fluctuations play a dual 
role in exciton transportjii£i2fi. In structures with large 
static disorder excitons are localized. Dynamic fluctu- 
ations in this case remove localization. However, if the 
fluctuations are strong this results in the dynamical local- 
ization of the exciton. A seminal model to treat these fast 
fluctuations is due to Haken, Strobl and Reineker^i ri22 ? 
where the dynamical fluctuations are described by real 



stochastic processes, e n (t) and V nm (t) and modeled using 
the density matrix^i or stochastic Schrodinger 124 equa- 
tions. In the original work the stochastic processes have 
been chosen as white noise, i.e. delta-correlated in time. 
This model has been shown to capture coherent and inco- 
herent excitation transfer on the same footing. Moreover, 
one should notice that the equation describing Haken- 
Strobel-Reineker (HRS) model for ID systems is equiva- 
lent to telegrapher's equation, which further reflects the 
similarity between the exciton propagation in molecular 
aggregates and a wave propagation in a medium. 

The Haken- Strobl- Reineker model, with several exten- 
sions, has been utilized to study exciton dynamics in nat- 
ural and artificial molecular aggregate o 120 i 123 ~ 22L. The 
drawback of this model is that the exciton population 
does not relax to thermal equilibrium. There have been 
many investigations and extensions which have included 
colored noise and tried to cure the thermalization is- 
sue, e.g j 128 i 129 . The original model also does not in- 
clude the effect of strongly coupled high frequency modes 
of the monomers which is essential to describe the vi- 
bronic structure present in many H-aggregates. The re- 
cent development of so called quantum state diffusion 
method a 13Q i 131 can be considered as an extension of the 
HSR model which cures these deficiencies at the cost of 
colored complex noise. 

Fig. fTTT a) illustrates how the HRS model can be ap- 
plied to exciton propagation on a brickstone lattice of 
TDBC dyes at ambient conditions 124 . The second mo- 
ment of the exciton wavefunction, Fig. fTTT b) shows that 
the initial propagation on the time scale of 20 — 30 fs 
is ballistic or wave-like, which is associated with the 
quadratic dependence of the second moment on time. For 
the longer timescale the time dependence is linear, which 
characterizes a diffusive motion. 



C. The vibronic model 

Let us first consider a single monomer, where only the 
two lowest electronic states are taken into account, see 
Fig. [21 In the Born-Oppenheimer approximation^ 2 - we 
write the nuclear Hamiltonian of the electronic ground 
and excited states as 



H% /e = lj2( P nj + U n /e (Q 



,0) 



(8) 



where Q n j are the (mass-scaled) nuclear coordinates of 
monomer n and P n j are the corresponding momenta. 
U 9 ^ e (Q n j) is the Born-Oppenheimer potential in the elec- 
tronic ground/excited state (see Fig. [2]). 

Then, for a molecular aggregate the Hamiltonian 
in the electronic ground state is given by H 9 = 

[J2n=i H nJ \9ei)(9e\\- We denote by V nm the transi- 
tion dipole-dipole interaction between monomer n and 
m, which, for simplicity, is taken to be independent of 
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FIG. 11: (a) Time dynamics of exciton wave function on a 2D lattice of TDBC molecules subject to static disorder (linewidth 
a — 70 meV) and dynamic fluctuations (linewidth r = 30 meV) of molecular electronic transitions. At time zero the exciton is 
localized in the center of the lattice. The contour plots show the population of lattice sites at four different times, (b) Second 
moments of the wave function corresponding to the exciton dynamics shown in (a). Two different transport regimes can be 
observed: the ballistic or a wave- like propagation on timescales of 20 — 30 fs (M^ ~ t 2 ) and the diffusive motion (M^ ~ t) 
at longer times. Reprinted with permission from Ref?^. Copyright 2012, American Institute of Physics. 



nuclear coordinates. The Hamiltonian in the one-exciton 
manifold, Eq. (jl]), is then given by 



TV 



TV 



H e = Y,W n \n)(n\+ J2 Vnm\n)(m\. (9) 



n=l 



n,m—l 



The potential surfaces are sketched in Fig. [2] for the case 
of one mode. This harmonic approximation is widely 
employed in the literature, e.gi&i&I 3 ^".! 3 ^ 



with the "collective" Born-Oppenheimer surfaces W n = 
H n + Em^n^m- Note that the structure of Eq. ® is 
very similar to the purely electronic Hamiltonian (J5j). 
The only difference is that now the energies e n are re- 
placed by operators for the nuclear motion. Upon trans- 
fer of excitation the nuclear wavepacket, according to 
Fig. [21 is no longer in its equilibrium position. Thus 
the excitation transfer is linked to excitation of nuclear 
motion. This vibrational model is sometimes combined 
with the static disorder approac h 111 ! 133 ! 134 . The vibra- 
tional model offers a clear idea regarding why the J-band 
does not exhibit vibrational structure and the H-band 
does. 

The coupling to vibrations usually slows down the 
propagation of the exciton, however, it can help to over- 
come energetic barriers caused by differences in the elec- 
tronic transition energies. In this way efficient directed 
transport along a biased chain can be achieved. 

Often it is justified to approximate the Born- 
Oppenheimer potentials, Eq. (|5J), as harmonic potentials 
where the surface of the electronically excited state is just 
shifted relative to that of the ground electronic state, i.e. 

we have H£ == \Y!jLi( P nj + u ljQlj) and for tne ex " 

e n + 2 Ej=l { P nj +^nj (Qnj ~ 

AQ n j) 2 ) where e n denotes the energy difference between 
the minima of the upper and lower potential energy sur- 
face and the vibrational frequency of mode j is denoted 
by uj n j and the shift between the minima of the excited 
and ground state harmonic potential is denoted by AQ 



2 2^j 

cited electronic state H\ 

\2 



D. Open system approaches 



Open system approaches are closely related to the vi- 
bronic model discussed in the previous subsection. In 
these approaches one typically considers an infinite set of 
harmonic oscillators forming the bath. It is easy to see 



that the Hamiltonian (J9j) can be written as H e = H, 



nj- 



sys T" 

# sy s-bath+#bath with H sys given by the purely electronic 
part (|5j) with transition energies e n that already include 
static overall shifts induced by the environment. The 
bath is taken as i^bath = J2 n ^-n an d the system-bath 
coupling is described by a linear coupling of the elec- 
tronic excitation of a monomer to the set of vibrational 
modes, i.e. # sy s-bath = ]T n | n n ) ( n m \ V. k nj Q nj . The 
coupling constant k n j describes the coupling of the ex- 
citation on monomer n to the harmonic oscillator with 
frequency uj n j. Typically, the frequency spectrum of 
bath oscillators is taken to be continuous so that the 
coupling constant becomes a continuous function of fre- 
quency. The open system approaches have been used 
to describe optical and transfer properties of light har- 
vesting system a 139 i 140 and can also be applied to explain 
properties of organic dye aggregates^. We note that 
the HRS model can be considered as a special case of the 
open system models for which the bath is Markovian (i.e. 
memoryless). 



E. Mixed QM/MM 

With the increase of computer capabilities, in recent 
years it has become possible to simulate aggregates start- 
ing from a microscopic description^ 4 ^ 2i^. Since a full 
quantum mechanical treatment is still out of reach, one 
uses mixed quantum classical approaches. In these ap- 
proaches usually the nuclei are propagated classically in 
the electronic ground state using molecular dynamics. 
Along these nuclear trajectories one calculates the (time- 
dependent) transition energies using quantum chemistry 
methods. These time-dependent transition energies can 
then be used to obtain input parameters for open system 
models or to use them directly as the energy fluctuations 
in HSR model. However, both approaches have their 
problems, since one is restricted to a classical propaga- 
tion in the electronic ground state. 
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FIG. 12: (a) Structure of a planar A/2 organic microcavity 
with a thin layer of J- aggregates embedded, (b) An example 
of a cavity photon and a J-aggregate exciton energy level anti- 
crossing. Reprinted with permission from Ref.^^. Copyright 
2010 by the American Physical Society. 



F. Semi-empirical approach based on 
polarizabilities 

While the approaches discussed above start from a mi- 
croscopic description of the exciton-phonon coupling and 
the environment, in the semi-empirical approaches one 
can take the standpoint that all the relevant information 
of the coupling of an electronic transition to other de- 
grees of freedom is already contained in the shape of the 
absorption spectrum. The absorption spectrum in turn 
is linked to the polarizabilities of the monomers. 

The induced electric moment of a monomer depends 
on the local field at that monomer, which is produced 
by the electric moments of all the other monomers and 
the external field. This leads to a set of coupled equa- 
tions from which one can extract, for a given arrange- 
ment of the monomers, the polarizability of the aggregate 
and thus the optical properties. Such equations have 
been derived using classical treatments^ or quantum 
approache s 117 ! 146 . This method implicitly assumes that 
the frequency dependent polarizability does not change 
(beside an overall shift in energy) when going from the 
monomer to the aggregate. Recently, it has been demon- 
strated that using measured monomer spectra as input 
one can accurately describe the band-shape of J- and H- 
aggregate o 60 i 147 . This approach also gives an intuitive 
explanation of the vibrational structure and the broad- 
enings of the absorption bands of the aggregate. 



ilar to the Auger effect, followed by a fast non-radiative 
energy conversion. Within an approximate picture this 
can be viewed as a process by which two excitons local- 
ized on neighboring molecules are combined into a higher- 
energy electronic excitation of a single molecule. Then, 
the resulting excitation quickly decays to the lowest ex- 
cited state or ground state dissipating the exceeding en- 
ergy through vibrations. A detailed kinematic model for 
the exciton-exciton annihilation in molecular crystals of 
different dimensions has been introduced by Suna ini 4 ^, 
and a similar master equation approach for exciton anni- 
hilation dynamics in natural photosynthetic complexes 
has been suggested in 149 . While the microscopic ap- 
proaches discussed above are still valid, the molecular 
Hamiltonian (J5j) should be extended to account for at 
least two excited states per each monome r 38 ! 15Q . 



V. HYBRID EXCITONIC STRUCTURES 

In this section we specifically focus on J-aggregates. 
Their large absorption and fluorescence cross-sections 
combined with a narrow line width and a small Stokes 
shift allow for a coherent coupling between excitons 
and photons in optical cavities or excitons and plasmon 
modes in metal structures. This property can be uti- 
lized as an interface between photonic, plasmonic and 
excitonic circuits. 



G. Beyond the single-exciton approximation 



A. Exciton-polariton structures 



Photonic device applications may require structures 
operating in a nonlinear response regime, where interac- 
tions between excitons cannot be neglected. For example, 
thin J-aggregate films switches based on optical bistabil- 
ity have been suggested in 3 -^. At higher intensities, how- 
ever, an additional loss-channel - exciton-exciton annihi- 
lation - appears. The underlying physical process is sim- 



Most of the studies of strong coupling between exci- 
tons in J-aggreagtes and photons have been done using 
organic micro cavities^i - -^. A schematic illustration of 
an organic microcavity - a planar A/2 cavity with a J- 
aggregate layer embedded in it - is shown in Fig. fl2T a). 
Strong coupling of cavity photons with excitons in a J- 
aggregate result in formation of exciton-polariton modes. 
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This is associated with the splitting of a cavity mode - 
vacuum Rabi splitting, see Fig. fT2l b). 

Within a simple model assuming only one exciton state 
and a single photon mode the value of the Rabi splitting, 
£1r, is just 



^r = dx • E Vi 



(10) 



where dx is the transition dipole of the narrow exciton 
transition, E vac is the vacuum field in the position where 
the J-aggregate is located, and we assume that the vac- 
uum field does not vary substantially on the scale of the 
aggregate. In J-aggregates the value of the exciton tran- 
sition dipole scales with the number TV* of coherently 
coupled molecules involved as dx ~ VN*. The value 
of N* is usually much smaller than N - the total num- 
ber of molecules in the aggregate, since the environment 
leads to localisation of the exciton states. The value 
of the Rabi splitting observed for J-aggregates strongly 
coupled with optical cavities ranges from several tens of 
me yi5i 1 i52 to severa i hundreds of me V 102 i 153 i 154 depend- 
ing on the disorder present in the aggregate and the de- 
sign of the cavity. Therefore, the splitting can be ob- 
served even at room temperatures. The photolumines- 
cence in these structures is usually observed from the 
lower polariton branch only, due to fast relaxation of po- 
laritons. By detunning the frequency of the cavity mode 
from the exciton resonance transition one could control 
the mixture of the exciton and photon and therefore mod- 
ify its coherence properties making it more photon-like 
or exciton-like. Electroluminescence from a polariton 
mode has been demonstrated using a light-emitting-diode 
structure with a J-aggregate layer^i Moreover, lasing 
from exciton-polariton mode has been shown for cavities 
with single anthracene crystals^ and some preliminary 
results were reported on lasing from cavities with cyanine 
dye J- aggregates^. Large interest has been attracted to 
the idea of polariton Bose-Einstein condensation (BEC) 
in organic cavities. This may bring the rich and con- 
troversial quantum physics of non-equilibrium polariton 
condensatesi££~"i££ observed at low temperatures in inor- 
ganic systemsiSiriSS up to room temperatures. However, 
to the best of our knowledge no confirmed observations 
of BEC in organic cavities have been reported yet. 



B. Plexciton structures 

A strong coupling of excitons in molecular aggregates 
with plasmons in noble metals has been demonstrated for 
both propagating plasmon modes in filmsi^ as well as 
localized modes in various types of nanostructure o 165 i 166 



and nanostructure arrays^. The hybrid plasmon-exciton 
modes were also named plexcitonsi^ 6 -. The observed val- 
ues of the splitting between the plexciton modes are of 
the order of several hundreds of meV. The linewidth 
of plasmon modes is usually sufficiently larger than the 
linewidth of excitons in molecular aggregates (plasmon 
lifetime is about tens of femtoseconds as compared to pi- 
cosecond exciton lifetime in J-aggregates). Therefore, the 
interaction between the exciton and the plasmon modes 
can frequently be considered as the coupling of a single 
mode (exciton) to a broader continuum (plasmon). This 
results in the formation of a Fano resonanc e) 168 ! 169 . It is 
interesting to notice that in the case of a strong optical 
pumping of plexcitonic structures non-linear Fano effects 
can be observed^ 



VI. CONCLUSIONS 

Aggregates of organic molecules - supramolecular as- 
semblies with strong resonant near-field interactions be- 
tween electronic transitions - could be exploited in the 
design of nanophotonic devices at the true nanometer 
scale. The molecules, forming the aggregates, interact 
collectively with optical fields, and the absorbed energy 
is transferred in the form of excitons on a sub-micron 
scale. The exciton transport within the aggregates pos- 
sesses coherent properties even at room temperatures, 
similar to the electromagnetic wave propagation through 
a medium. Moreover, molecular aggregates can be cou- 
pled coherently to photonic and plasmonic structures. 
While at the present stage there is a sufficient gap be- 
tween the research communities studying excitonics and 
photonics, this review calls for merging the knowledge 
from the two fields. 
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